Resonant origin for density fluctuations deep within the Sun: 
hehoseismology and magneto-gravity waves 

CP. Burgess^rN.S. Dzhalilov^rx.I. Rashba^'^* V.B. Semikoz^'^* J.W.F. Valle^* 



cn 

O 

o 

< 

00 
(N ■ 

> 

^ ■ 

o 
m 
o 

Oh 

6 



^Physics Department, McGill University, 3600 University Street, Montreal, Quebec, Canada, H3A 2T8 
^ The Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation of the Russian Academy of Sciences, 

IZMIRAN, Troitsk, Moscow region, 142190, Russia 
^AHEP Group, Instituto de Fi'sica Corpuscular, C.S.I.C./Universitat de Valencia 
Edificio de Institutos de Paterna, Apartado 22085, E^6071 Valencia, Spain 



IFIC-03-12; McGill-02/40 

Key words: Magnetohydrodynamics (MHD) - Sun: helio- 
seismology - Sun: magnetic fields - Sun: interior - neutri- 
nos 



Abstract 

We analyze helioseismic waves near the solar equator in 
the presence of magnetic fields deep within the solar ra- 
diative zone. We find that reasonable magnetic fields can 
significantly alter the shapes of the wave profiles for he- 
lioseismic g-modes. They can do so because the exis- 
tence of density gradients allows g-modes to resonantly 
excite Alfven waves, causing mode energy to be funnelled 
along magnetic field lines, away from the solar equatorial 
plane. The resulting wave forms show comparatively sharp 
spikes in the density profile at radii where these resonances 
take place. We estimate how big these waves might be in 
the Sun, and perform a first search for observable conse- 
quences. We find the density excursions at the resonances 
to be too narrow to be ruled out by present-day analyses 
of p-wave helioseismic spectra, even if their amplitudes 
were to be larger than a few percent. (In contrast it has 
been shown in 1 8 1 that such density excursions could affect 
solar neutrino fluxes in an important way.) Because solar 
p-waves are not strongly influenced by radiative-zone mag- 
netic fields, standard analyses of helioseismic data should 
not be significantly altered. The influence of the mag- 
netic field on the g-mode frequency spectrum could be used 
to probe sufficiently large radiative-zone magnetic fields 
should solar y-modes ever be definitively observed. Our 
results would have stronger implications if overstable solar 
g-modes should prove to have very large amplitudes, as has 
sometimes been argued. 
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1 Introduction 

Helioseismology has become a precision tool for studying 
the inner workings of the Sun, providing one of the only 
direct probes of physical properties as a function of depth, 
below the convective layer and into the radiative zone. The 
frequencies of thousands of normal modes have been mea- 
sured, and current levels of agreement between these mea- 
surements and calculations preclude deviations of density 
profiles and sound speeds (as functions of depth) in ex- 
cess of one percent from solar-model predictions. Careful 
analysis of small discrepancies between models and mea- 
surements which have arisen in the past has helped refine 
the models, by winnowing out small errors in opacities and 
other quantities used as inputs. 

The vast majority of heUoseismic analyses to date are 
performed in the approximation which neglects the sun's 
magnetic field. This approximation is generally very rea- 
sonable since the energy density, B^/Stt, of the expected 
fields is much smaller than other energies in the problem, 
such as gas pressures, deep in the solar interior This is 
expected to be particularly true in the solar radiative zone, 
below the turbulent convection which is believed to be re- 
sponsible for the solar dynamo. 

Very little is directly known about magnetic field 
strengths within the radiative zone - see 1131 [D for early 
studies. A generally-applicable bound is due to Chan- 
drasekhar, and states that the magnetic field energy must 
be less than the gravitational binding energy: B'^/Stt < 
GMl^/R^, 01 B < 10^ G. Sti-onger bounds are possible 
if one makes assumptions about the nature and origins of 
the solar magnetic field. For instance, if it is a relic of the 
primordial field of the collapsing gas cloud from which the 
sun formed |27 1, then it has been argued that central fields 
cannot exceed around 30 G |6|. Similarly, the dynamo 
mechanism can only generate a global field in the radiative 
zone of a newly -born Sun with amplitude below 1 G II21I . 
Even stronger limits, B < 10~^ G apply |25| if the solar 
core is rapidly rotating, as is sometimes proposed. On the 
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other hand, it has recently been argued 1151 that fields up to 
7 MGauss could persist in the radiative zone for billions of 
years and are consistent with current observational bounds. 
Other authors 1121 have recently entertained radiative-zone 
fields as large as 30 MGauss. Since the initial origin and 
current nature of the central magnetic field is unclear, we 
consider as admissible any magnetic field smaller than of 
order 10 MGauss. 

The purpose of the paper you are now reading is to ex- 
amine the influence of magnetic fields in more detail, deep 
within the solar radiative zone. Our motivation for so do- 
ing is threefold. First, although magnetic fields in the ra- 
diative zone are believed to be small, they cannot be di- 
rectly measured, and so quantitative bounds on their size 
require the explicit calculation of their effects. Ours is a 
first step along these lines, where we assume a particularly 
simple geometry which we argue to be appropriate when 
the background magnetic field is orthogonal to the local 
density gradients. Such configuration would occur near 
the equatorial plane in a dipole magnetic field. Second, 
plasmas in magnetic fields are notoriously unstable, with 
small perturbations often giving rise to strong effects, per- 
haps amplifying the implications of a small magnetic field. 
Indeed our analysis uncovers such instabilities, which give 
rise to surprisingly large effects. Third, the recent advent 
of neutrino telescopes opens a new window on the solar 
interior, motivating the exploration of the implications of 
solar magnetic fields for solar-neutrino properties. 

In a nutshell, our main findings are these: 

• We find, contrary to naive expectations, that pro- 
files of helioseismic g-modes as a function of solar 
depth can be appreciably altered in the radiative zone 
even for magnetic fields which are not unreasonably 
large 1271 l6l 1211 . In particular, density profiles due 
to these waves tend to form comparatively narrow 
spikes at specific radii within the sun, correspond- 
ing to radii where the frequencies of magnetic Alfven 
modes cross those of buoyancy-driven (g-type) helio- 
seismic modes. Due to this resonance, energy initially 
in g-modes (presumably due to turbulence at the bot- 
tom of the convective zone) is directly pumped into 
the Alfven waves, causing an amplification of the den- 
sity profiles in the vicinity of the resonant radius. This 
amplification continues until it is balanced by dissipa- 
tion, resulting in an unexpectedly large density excur- 
sion at the resonant radii. 

• Since the only such level crossing which occurs 
within the radiative zone occurs with g-modes, and 
not p-modes, the latter are less substantially affected 
by radiative-zone magnetic fields. This makes it very 
unlikely that radiative-zone magnetic fields can sig- 
nificantly alter standard analyses of helioseismic data, 
which ignore solar magnetic fields. 

• We identify the implications of magnetic fields for the 
amplitudes and frequencies of solar g-modes, in prin- 



ciple permitting inferences about radiative-zone fields 
to be drawn once these modes are experimentally ob- 
served. We find that the bound implied by the consis- 
tency of our analysis with existing discussions of p- 
mode heUoseismic waves is at present uninterestingly 
weak. 

The main message in these findings is that the detailed 
shape of helioseismic g- waves can be significantly changed 
by reasonable radiative-zone magnetic fields, and this sen- 
sitivity to magnetic fields is driven by the occurrence of 
level crossing between helioseismic g-modes and Alfven 
waves. Of course, although this cartoon describes the 
Alfven and helioseismic waves as if they were separate, 
a real calculation should simply diagonalize the normal 
modes of the linearized hydrodynamic problem including 
the magnetic fields. We perform such a diagonalization 
explicitly, and find the magneto-gravity (MG) waveforms 
whose peaks correspond to the above resonance argument. 
This diagonalization allows us to determine how the mode 
energies and excitation rates depend on the strength of the 
background magnetic field. 

The resonance we identify has also been recognized as 
possibly playing a role elsewhere in the Sun. In particu- 
lar, the feeding of energy into this resonance has been pro- 
posed I iF 19" '2(7, 16", 171 as a mechanism for heating the 
solar corona. Indeed the analysis of these authors is very 
similar to the one we present here, with the main difference 
being their neglect of gravitational forces (as is appropriate 
for applications to the solar atmosphere, but not for the ra- 
diative zone). 

Indeed, we believe the study of the implications for neu- 
trinos of the Alfven resonant layers is well worth pursu- 
ing, particularly since density fluctuations themselves can 
have effects, without invoking the existence of large mag- 
netic moments. Two other lines of argument also encour- 
age such a study. First, previous experience along these 
lines 1 22 23 29 , 28 , 24 26 T| argues that neutrinos are 
most sensitive to density fluctuations near the depths at 
which MSW resonant oscillations take place, r ~ 0.3 Rq, 
where the resonances we derive indeed can arise. Second, 
the density spikes we here find are much narrower than 
are generic helioseismic g-wave profiles, possibly permit- 
ting Alfven waves to cause effects which garden-variety 

waves cannot [4J. We present a discussion of neutrino 
propagation through these waves in Ref. (8), where we find 
potentially interesting prospects for detectable effects if 
radiative-zone magnetic fields lie in the range 10-100 kG. 

We organize our presentation as follows. In section 2 
we set up the hydrodynamic problem in the presence of 
gravity and magnetic fields, and give the total set of MG 
equations whose solution is required. We next derive the 
linearized MHD equations, which we solve by relating all 
physical quantities to the radial component of the perturba- 
tion magnetic field amplitude. In section 3, we apply rea- 
sonable boundary conditions to set up an eigenvalue prob- 
lem which determines the dispersion relation for the dis- 
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Crete eigenfrequency spectrum, a;„(i?o)- For a simple ge- 
ometry we identify the labels, n, and spectrum of the MG 
eigenwaves. In section 4 we focus on the resonance layer 
position and compute how its position and width vary from 
mode to mode. We also here present results for the size 
of the maximum density excursions at the resonant points 
as a function of magnetic field strength. Finally, section 5 
discusses some observational constraints on the wave fre- 
quencies and amplitudes. 

2 The master equation for MG 

waves 



V^(/) — —AnGp. As usual, G here denotes Newton's con- 
stant. The last term of equation (|3} expresses the contribu- 
tion of the Lorentz force to the local momentum budget. 
Finally, the system is completed by Faraday's equation. 



— = rot [v X B] 



(4) 



that governs the time evolution of the magnetic field. Here 
Ccond represents the fluid's conductivity. In what follows 
we take the plasma to be an ideal conductor, meaning we 
take acond to be large enough to neglect the last term in 
this last equation. Clearly these equations reduce to those 
of standard Helioseismology Models (SHM) in the hmit 



In this section we briefly review how magnetic fields in- 
fluence the equations of hydrostatic equilibrium on which 
helioseimic analyses are based. To this end we first set up 
the relevant MHD equations, and linearize them to iden- 
tify the differential operator whose spectrum dictates the 
allowed MG wave frequencies. 

Two equations of fluid dynamics are not changed by the 
presence of a magnetic field - at least within the regime of 
present interest. The first of these is the continuity equa- 
tion, expressing conservation of mass: 



dp 
'dt 



+ pu — , 



(1) 



where d/dt = d/dt + v • V is the usual convective deriva- 
tive, with V representing the fluid velocity, and p and p 
are the fluid's mass density and pressure. The variable 
u — div V need not vanish if the fluid is compressible ^ 

The second unchanged equation expresses energy con- 
servation for a polytropic medium. This equation states: 



dp 
'dt 



pdp 
p dt 



-(7-l)Q 



(2) 



where 7 — Cp/cy is given by the ratio of heat capacities, 
and Q is the sum of all energy density sources and losses. 
In principle, electromagnetic processes enter into Eq. (|2jl 
through their contributions to Q, but for ideal MHD we 
neglect both the heat conductivity and viscosity contribu- 
tions to energy losses, as well as the ohmic dissipation, 
Q = j'^ /<^cond, where acond is the electrical conductivity. 

The magnetic fields enter more directly into the Eu- 
ler equation (conservation of momentum), which is of the 
form 



dv 1 



rot B X B 



(3) 



with the local force of gravity given by pg. The local ac- 
celeration due to gravity is related to the Newtonian po- 
tential by g = V0, with (p given by the Poisson equation 



' This variable u was the main focus of our earlier asymptotic analysis 
of MG waves 14J . 



2.1 Linearized form 

Our interest is in small-amplitude oscillations about a 
background configuration, and so at this point we split 
all variables into background and fluctuating quantities, 
A = Aq+A', with the A' denoting small fluctuations about 
the background value Aq. Since we take our background 
configuration to be static, vq — 0, we regard v to be pure 
fluctuation, but suppress the prime on v in the interests of 
notational simplicity. 

The system of linearized equations obtained in this way 

is 



9p' , 
at 



0, 



9v 1 r 

Pott- + Vp' - gp' - — rot B' x Bo 
at in I 

dp' 

+ (v • V)po + lPou = , 

dB' 
~dt 



= 0, 



rot ( V X B, 



») 



(5) 



where, as described above, v, p', p' and B' are the small 
Euler perturbations. 

In addition to our previously-mentioned assumptions of 
ideal conductance - i.e. Ccond 00 - and static back- 
grounds - Vq = 0, with background quantities time in- 
dependent - these equations make three further assump- 
tions. First, they adopt the Cowling approximation, which 
amounts to the neglect of perturbations of the gravitational 
potential, (i.e.: (f)' = 0). Second, they assume the fluctua- 
tions are adiabatic, Q' = 0, as is satisfied with good accu- 
racy by high frequency oscillations. Finally, they assume 
the background field to satisfy rot Bq = 0. 

The backgrounds poiPo etc., must themselves satisfy 
the hydrodynamic equations, which for negligible (or con- 
stant) magnetic fields reduce to the usual equations of hy- 
drostatic equilibrium, Vpo = Pog- Finally, we suppose 
there to be no energy sources or sinks in the background, 
Qo = 0. Under these circumstances the background den- 
sity po is as given by Standard Solar Models (SSM) (2||3l. 

If we also assume 7 to be constant - as is a good approx- 
imation within the whole radiative zone - we can derive 
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from the system of linearized equations (|5} an evolution 
equation for the velocity field v: 



= V(g-v)+(7-l) gu+c^Vu-t-^ _ 



— X rot b 
Bo 

(6) 

1/2 

where b = B'/Sq and va = Bq/ (47r/9o) defines the 
Alfven speed. 

2.2 Geometry 

To proceed further we next choose the geometry of the so- 
lution we shall seek. Our interest is in a magnetic field 
which is perpendicular to the local density gradients and 
gravitational fields. Our interest is also in the deep interior 
of the radiative zone, but not right down to the solar center. 

For these purposes we can take Bq to be both constant 
and approximately uniform. It then suffices to work within 
an approximately rectangular, rather than cylindrical, ge- 
ometry. In these circumstances it is convenient to choose 
a Cartesian coordinate system whose z-axis is the 'radial' 
direction, as defined by (but opposite in direction to) the lo- 
cal gravitational acceleration, g. With this choice we take 
z e (0, Rq), where z = represents the solar center and 
z = Rq denotes the solar surface. We are actually mainly 
interested in the radiative zone, for which z < 0.7 Rq. 

Since we have in mind magnetic fields perpendicular to 
physical gradients and gravitational fields, we take Bq to 
lie along the x-axis and the gradients to lie along the z- 
axis, leading to the components Bo = {Bq, 0, 0), and g = 
{0,0, ~g{z)). 

2.3 Isolating variables 

Returning to our original problem, we wish to solve the lin- 
ear MHD equations without resorting to the WKB approx- 
imation. Our first goal to this end is to eliminate all vari- 
ables in terms of a single one, which we choose to be the 
z-component of the perturbed magnetic field, bz{z). The 
spectrum of MHD modes is then found from the eigenval- 
ues of the linear ordinary differential equation satisfied by 

This is most easily done by transforming v and b to 
Fourier space in the x and y directions, such as in v = 

v(z) exp^— i(a;t — k^x — kyy)j . In this way we find that 
the velocity perturbations are given by the system 









Vy 


1 


—ik 




1 




Vz 




oz 



ikx{clu - gvz) 



-iky{(?gU — gVz) — iuiv\ rot ^ b 
) + (7-l)5u-c2|^ + 



(7) 

, (8) 

(9) 
(10) 



+ i'jjv\ rot j,b] . 
The equation for the magnetic field perturbation b(z) is 



then derived from the last equation of eqs. (|5}, 

iwb = ikx^ — ^ u , (11) 

Bq 

and by the definition of the compressibility, u: 

dv 

u = ik^Vx + ikyVy + . (12) 

In our later applications we are most interested in the 
density perturbation given in terms of the other quantities 
by 

hz+u. (13) 



. P 

Po k^H 



Substituting the expressions for Vx,y from eqs. (0 and 
(I12> . and finding the magnetic field component by from the 
Maxwell equation div b = we find expressions for the 
dynamical variables in terms of bz. In so doing we also 
use two further approximations. We first assume the adi- 
abatic parameter, 7, to be constant, and we also use the 
low-frequency approximation, ^ k^cl. (This latter 
assumption has the effect of filtering out the acoustic p- 
modes from our analysis.) We find in this way: 



bx{z) 

by{z) 

Vz{z) 
Vx{z) 

Vy{z) 

Po 
Po 

where k'j^ - 
becomes 

u{z) = 



ikx dbzjz) 
k^ dz 
iky dbz{z) 
dz 

hziz) , 



kl 



k_ 
ikxUj 



ikyijj 
1^ 



70; 



gbzjz) 

kxC^ 

gbzjz) 

kxcl 
dbz{z) 



u? \ 1 dbz[z) 



kid 



1,2 2 



1 dbz{z) 



kj- dz 



i+i) /7-1 



bz{z) + 



uj^H dbz(z) 
k j^c^ dz 



kxH \ 7 

kl + ky. Using these, expression (fT2l for u 



(14) 



OJ 

kx 



g 1 dbz{z)\ 

-bz{z) - ——J— ) ^0 



dz 



(15) 



Note that we retain here in the transverse velocity compo- 
nents, Vx,y, terms which are subdominant in the quantity 
w^/fc^c^ <C 1, since these are needed to obey the com- 
pressibility condition, Eq. J12> or (I15> . 



3 The eigenvalue problem 

We may now identify the equation satisfied by bz{z) it- 
self, whose solution determines all of the other variables 
through the expressions derived above. We find bz to be 
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determined as the solution to the following second order 
linear ordinary differential equation, 

/ klvl\ d^M N'db.jz) ^ 
\ Lo"^ ) dz'^ g dz 

+kl('^-l + ^)bAz)^0, (16) 

where N denotes the Brunt-Vaisala frequency, as defined 
by 

N (z) = g{z) — . (17) 

V7P0 dz po dz J 

Equation ( I16t is one of our main results, and the remainder 
of the paper is devoted to the construction of its solutions. 

3.1 Qualitative Analysis 

Eq. ( I16> describes the interaction of adiabatic gravity 
modes with magnetic fields (magneto-gravity modes), in 
the approximation of low frequencies, oj^ ^ k'^c^- Be- 
fore pressing on, it is instructive to examine the qualitative 
properties of the solutions to Eq. ilbt . since these capture 
the results we will obtain from a more detailed analysis. 

In the limit Bq ^ (and so va 0), Eq. ( I16> re- 
duces to the standard evolution equation for 'pure' helio- 
seismic g-modes, which is usually expressed in terms of 
the variable (which is related to bz through the relation 
Vz = uibz/kx). Since more detailed analyses - see Fig.^- 
show that TV rises from zero at the solar center, remains ap- 
proximately constant N « iVo, through the radiative zone, 
and then falls to zero again at the bottom of the convective 
zone, these g-modes can be thought of as the eigen-modes 
of oscillations inside the cavity in between the two regions 
where TV goes through zero. For a given wave frequency, 
u! < Nq, the turning points of this cavity are given by the 
condition oj ^ N. For smaller uj the lower turning point 
gets closer to the solar center, z ^ 0, and the upper one 
gets slightly closer to the bottom of the convective zone 
(CZ). 

Conversely, if gravity is turned off (N 0), then 
Eq. il6\ describes Alfven waves, which oscillate with fre- 
quency u! = kxVA and propagate along the magnetic field 
lines. Notice that since va oc p^^^^ this frequency grows 
with z, since the density of the medium falls. 

Keeping both magnetic and gravitational fields intro- 
duces qualitatively new behavior, as may be seen mathe- 
matically because equation il6\ acquires a new singular 
point which occurs when the coefficient of the second- 
derivative term vanishes. Since this singularity appears at 
the Alfven frequency, 

LO = t^VA (18) 

it can be interpreted as being due to a resonance between 
the g-modes and Alfven waves. This resonance turns out to 
occur at a particular radius because the Alfven frequency 
varies with radius while the g-mode frequency is indepen- 
dent of radius (see Fig.Q. The resonance occurs where the 



growing Alfven mode frequency crosses the frequency of 
one of the (/-modes, and the resulting waveforms would be 
expected to vary strongly at these points. 

The resonance can occur inside the radiative zone if the 
Alfven frequency climbs high enough to cross a g-mode 
frequency before reaching the top of the radiative zone. 
Since va oc Bq, whether this occurs or not depends on the 
central field value, Bq. Our main task is to find the depen- 
dence of the spectrum of eigen-modes on magnetic field, 
and to elucidate when such singular resonances can appear 
within the Sun. Of particular interest for observational pur- 
poses is to know how the density and sound speed vary 
near the resonances, especially when the resonance region 
occurs near the locations of resonant neutrino oscillations. 

The presence of a resonance at a particular radius, z = 
Zr, shifts the position of the upper turning points of the cor- 
responding magneto-gravity wave towards deeper regions 
of the Sun. This leads to several new effects. First, the 
shortening of the cavity due to the presence of magnetic 
field causes the eigen-frequencies to depend on the mag- 
netic field value. Second, there can be energy transfer be- 
tween g-modes and Alfven waves, within the narrow sin- 
gular resonance layer, leading to the corresponding eigen- 
frequencies acquiring imaginary parts. Third, the nearer 
the upper boundary of MHD cavity is to the solar center, 
the stronger the g-modes are confined to the solar core. 

More details of these waveforms can be obtained from 
a WKB-type analysis of the master Eq. il6\ . Near the so- 
lar center where iV^ 0, if {kxVA)'^ / oj'^ <C 1 one ob- 
tains the exponential solution bz e^^^^ . (Which com- 
bination of these solutions appears is fixed by boundary 
conditions, e.g. if 6z(0) = at the solar center we have 
bz oc sinh(A:^z).) One finds similar behaviour for the so- 
lution above the singular resonant layer, z > z^, up to the 
top of the radiative zone. This exponential growth happens 
because of the exponential decrease of density (and so ex- 
ponential increase in va) with z. The requirement for com- 
plex frequencies arises from the demand that the solution 
remain regular in the narrow Alfven resonance layer. 

3.2 Background Configurations 

In order to acquire analytic solutions to eq. ( II 6t . we must 
specify a background profile for pq{z) and g{z). We model 
Poiz) as an exponential density profile, po{z) = Pc e~^^^ 
with H constant, and take g to be approximately constant. 
We now argue these to be reasonable approximations well 
inside the radiative zone, and not too close to the solar cen- 
ter 

The variation of background quantities may be charac- 
terized by the scale heights, HJ^^ := AQ){dAQl dz)\, 
for A = p.,p,T etc.. Solar models G1|3|2Q| show that the 
density scale height is the shortest in the radiative zone, 
Hg, Ht > Hp, and that Hp takes the constant value 
Hp ^ H ^ 0.1 Rq to a good approximation, except near 
the solar center and solar surface. On these grounds we 
have that the gravitational equation, dg/dz = —AirGpo, 
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Figure 1 : Relevant frequencies plotted against radial position within the sun. The solid curve gives the Brunt- Vaisala frequency, 
N{z), while the dot-dashed lines are the Alfven frequencies. Notice that N{z) goes to zero at the solar center and the top of 
the radiative zone. The horizontal dashed lines represent the frequencies of magneto-gravity waves computed from Eq. i35\ 
for several choices for the mode number, n, for fixed values Kx — K± = 100. As discussed in the text, resonances occur 
where the mode frequencies intersect the Alfven frequency, indicated by circles in the figure. The two panels correspond to 
two choices for the magnetic field: 10 kG (left) and 100 kG (right). Minimum and maximum mode numbers are indicated, with 
riniin defined by the condition that its resonance occurs at the top of the radiative zone, Zr{nynin) — 0.7Rq, and n,nax having 
resonance at Zrl^^max) = O.12i?0 (see Figs.|4]and|5}- The dotted Unes denote acoustic (Lamb) frequencies for I — 10 and 100 
(see 1 1 1 1) while the double-dot-dashed horizontal line represents the trapping region for a p-wave of frequency 3000 /iHz and 
I — 10. The horizontal solid line represents the trapping region for a g-wave of frequency 100 /iHz. 
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implies g — AirGpcH (l — e^^/^), which is approxi- 
mately constant for z much larger than H. For constant 
g the hydrostatic equation, dpo/dz — —po g, then implies 
Po = PcC^p{—z/H) which also implies constancy of the 
temperature for gases obeying the ideal gas equation state, 
Po/ Po = RTq/p- Here R is the gas constant, R = kg/mu, 
fjL is the molecular weight measured in atomic mass units 
m„ ~ mp, and kg is Boltzmann's constant. In terms of 
these quantities the density scale height becomes: H — 
(P-jg (= constant) d 13, with c = ^/KI\JJl being the 
isothermal sound velocity, which is related to the adiabatic 
sound velocity, c^, by = . All of these quantities are 
constants within the approximations we are using. 

Notice that our approximations that j, H = O.li?© and 
g are all constant also imply approximate constancy for the 
Brunt-Vaisala frequency, N, within the radiative zone. 

3.3 Solution Using Hypergeometric Func- 
tions 

Changing variables to ^ = k'^v\/ijj'^ ~ e^/^ we recast 
Eq. M6\ as 



2 d%.{^) 



de 



7-l^ db^io 



1 



(19) 



where we introduce the dimensionless parameters K'j^ = 

Since the equation ( I19> has three regular singular points 
- at ^ = 0, 1 and cx3 - it can be put into standard hyperge- 
ometric form. 



dY 

m — -ahY = ^, (20) 
di 



through a change of the dependent variable = Y(f^- 
This permits an explicit solution for 6, in terms of Gaussian 
hypergeometric functions, F(a, 6; c; = 2-^1(0, 5; c; 0. 
Here the hypergeometric coefficients, a, h and c, are given 
by 

a = a + K^, b = cT-Kj_, c^2a + j^^, (21) 
where the complex index, cr, is defined by the expression 

7-1 



with 



2a = 



7-1 



7 



q , 



(22) 



- AKKk^ - 1) 



1/2 



General Solutions 



Two linearly-independent solutions to the Hypergeometric 
equation, (I20> . are 

Yi = F{a, b; c; , ^2 = ^^-^F{a~c+1, 6-c+l; 2-c; C) . 



These solutions must be analytic for all ^ ^ 0, 1 and 00, 
but they may acquire branch points at these three regular 
singular points depending on the values of a, b and c. In 
the present case because Rc (c — a — 6) = 7^^ > 0, the 
solutions are bounded at ^ = 1 as long as neither c nor 
2 — c is zero or a negative integer. Moreover, the second 
fundamental solution, Y2, is finite at ^ = if Rc (1 — c) = 
Re g > 0. 

We shall find that the waves of interest have complex 
frequencies, which we parameterize as = llJi{1 + id) 
with LUi and d real, with \d\ <^ 1. Complex uj implies, 
from i22\ . that the parameters q (and so also c) are also 
complex, q ^ qi + iq2, with real and imaginary parts: 



q = qi + iq2 = ^/wi + iw2 = 



1 

7i 



\j \J w\ + ~u^ + wi + i sgn(d) ]J ^wj + ~w^ - wi 

(23) 



where wi^2 are given by: 

2 



Wl 



7-1 
7 



(1 - 



ufii + d-^y 



W2 



SdN^Kl 



tj2 (1 + ^2)2 • 

It is clear from these equations that one always has Rc q = 
91 >0. 

The general solution for bz{z) is therefore given by 
the linear combination of the independent solutions to 
Eq. (|20j: 

b,iz)^D,CFia,b;c;0 + 

+D2 ^i-'^-i/7^(i + a _ c, 1 + 6 - c; 2 - c; ^^4) 

where one combination of the integration constants, 
D1/D2 is obtained by imposing boundary conditions at 
= ^{z = 0) and ^* = ^(z = z*), where « O.7i?0 
denotes the top of the radiative zone. 

Resonance Conditions 

Although the solutions are bounded at the singular point 
^ = 1, they are generally not analytic there. The position 
Zr determined by the condition 



Re^(zr) = {kxVA{zr) /cuif 



1 



(25) 



therefore defines the position of the 5-mode/Alfven reso- 
nance. This resonance occurs within the radiative zone 
only if < 1 and > 1. In this section we ask what 
this requires for the properties of the wave. 

Expressed in terms of the central quantities, the Alfven 
resonance condition is Re ^z^) = 1, or Re = 
(k^VAc/t^i)^ = e""^--/^. Since < z^ < w O.Ti?© 
this condition limits the values of K± for which resonances 
can occur. Using the numerical values _ff = 7 x 10^ cm, 
VAc = (-B0/43.4 G) cm s^^ (obtained using the central 
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density Pc = 150 g cm^''), N = 10^ (obtained using 
the mean acceleration g = 0.2 km s~^ and the ideal-gas 
adiabatic parameter 7 = 5/3), we find the following limits 
on the transverse wave number Kji_ > ~ kxH: 



7 X lOn^i/iV) 3P 
- if,(i?o/43.4G) 



(26) 



where the lower (upper) limit corresponds to the choice 

Zr = (Zr = ~ 0.7i?o). 

For frequencies, uJi ~ N, and moderate central mag- 
netic fields, Bq < 100 G, we see that the existence of a res- 
onance requires the wave number K± to be large, K± ^ 
1. Smaller K± necessarily requires either larger mag- 
netic fields or extremely low frequencies or both. Since 
we cannot, within the scope of our approximations, re- 
duce the angular frequency lower than the solar angular 
rotation frequency (for which the Solar rotation period is 
Tq ~ 27 days), a reasonable lower limit for uji might be 
(jji ~ 10~^ s~^ (corresponding to a period 2tt/uji ~ 8 
days). Even for such small frequencies one finds from 
Eq. (|26} the condition 3 x lO'^ > i^x > 10^, as required. 

Boundary Conditions 

At the solar center (z = and so ^ = ^c) we use the 
boundary condition which would have arisen as a smooth- 
ness condition in cylindrical coordinates: 



This requires the integration constants to satisfy 



£1 ^ 

Di ^'^ F{1 



F{a,b; c;^c 



c, l + fe-c;2-c;^c) 



(27) 



(28) 



This boundary condition is sometimes called the 'total 
reflection' condition because of the observation that it im- 
plies (for real uj) that the reflection coefficient at z = 
equals unity: R = \D2/Di\'^ = 1. To see why this is so 
notice that real uj implies q is pure real or pure imaginary. 
Within the domain of the approximation k'^cl/uj'^ ^ 1 it 
is pure imaginary, so we write q ^ —i(3 for some real (3. 
With this choice we have cr — > (7 — l)/(27) + i/3/2, and 
so 

a-^a + K^^ b ^ (T — K±, and c — > 1 + i/3 , 
which leads to the relations: 

1 + a — c 



l + b-c 
2-c 



27 2 

l-iP ^ c*. 



(29) 



From these conditions follows the conclusion R ~ 1. 

We next turn to the boundary condition at the top of the 
radiative zone, « 0.7 Rq. Since our interest is in res- 
onant waves which lie deep in the radiative zone, we wiU 



assume ^» = ^(2^*) S> 1 when applying this boundary con 
dition. 

The behaviour of bz for ^ > 1 is found using Eq. (|62] 
from Appendix A together with Eq. (I24> : 



xF ^a,l 

A2 - i-TAi 



c + a;l + a 
D2 



1 



Di 



xF 1 



b,b;l 



, 1 



where 



Ai 
A3 
A4 



nc)T{b-a) ^ 
T{b)T{c-a)' ' 

r(2- c)r(6- a) 
r(6-c+l)r(l-a)' 

r(2-c)r(a-6) 

r(a- c + i)r(i - b) ' 



r(c)r(a - b) 
r(a)r(c-6)' 



(30) 



(31) 



For ^ ^ 1, the two terms in Eq. i30\ behave as ^^^^ oc 
g±K± z/H ^ ^jj^ either fall or grow exponentially as func- 
tions of z. We take the absence of the growing behaviour 
to be our boundary condition for ^ = This requires 
vanishing of the whole factor within brackets in the sec- 
ond line of Eq. ( I30> . thereby providing a second condition 
which must be satisfied by the ratio D2/ Di. 



£2 



r(c)r(i-6)r(i 



c) 



r(2-c)r(a)r(c-6) 



(32) 



Given this boundary condition, b^ falls exponentially for 
z > Zr, bz ~ = S.Q^^^e^^^^^", as given by the 

first fine in Ea.( l30> . 



3.4 Eigenfrequencies 

The consistency of the two conditions, eqs. J32t and J28> . 
for D2 /Di gives the eigenvalue condition which fixes the 
dispersion relation for the resonant waves: 



F{a,b; c;^,) 



F{1 + a - c,l + b ~ c;2 - c;^c) 

^ r(c)r(i-b)r(i + a-c) 

^ ' r(2-c)r(a)r(c-5) ' ^ ' 

from which we read off the MHD frequency, 
u!{kx,ky,Bo,n), where the integer n is the mode 
number for the eigenfunction of interest. The solutions 
which are obtained in this way have the important property 
that they are complex, cu = wi (1 + id), a feature which 
is also present in the absence of gravitational loading, and 
which has been exploited in other contexts in attempts to 
explain the nature of coronal heating 1.18. .19. .20.. 16.. 17J . 
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We note in passing that, in contrast to the master equa- 
tion in ( I16> . Eq. ( I33> does not reduce to the standard g- 
mode spectrum in the Hmit Bq — 0. It does not do so 
because of our assumption that N be approximately con- 
stant. Although this approximation suffices for the study 
of MG waves, it is too crude to capture the resonant fre- 
quencies of g-modes. It cannot do so because g-modes are 
trapped inside the radiative zone precisely because N is not 
constant, but instead drops to zero at the top of the radia- 
tive zone. For MG waves this dropping of N is not crucial 
because its role in trapping the modes is instead played by 
the resonant layer at z — z^. 

The eigenvalue condition can be written in a more trans- 
parent way by taking advantage of some properties of hy- 
pergeometric functions. Recall that our interest is in the 
regime K±N/u!i ^ 1, and so the parameter q, as de- 
fined by Eq. i22i . is mostly imaginary, q = —i/3, since 
|d| < 1, with (3 = 2Kj_N/u; > 1. This implies that 
cr = (7— l)/(27)+i/?/2 is also large inasmuch as |o-| ::§> 1, 
and so must also be a and b. This observation allows us to 
use the Watson asymptotic form for hypergeometric func- 
tions (I65> - see Appendix A - which applies for large val- 
ues of a and b. 

Using this, together with liberal use of the convolution 
identity r(2;)r(l — z) = tt/ sin(7rz), allows us to rewrite 
the eigenvalue condition, Eq. (I33> . as a simple transcen- 
dental equation: 



2i/31n?7 = 2iTm + In 



tan ■ 



7 



± in 



1-7 
7 



(34) 



where the upper (or lower) sign on the right-hand- side 
is chosen if the imaginary part of the frequency is posi- 
tive: d > (or negative, d < 0). The variable 77 is de- 
fined hy r] — ^c^^^ — — 1, and depends on the 
magnetic field Bq through the Alfven velocity vac, since 
£,c — {kxVAc/^iY ■ The integer n — 0, ±1, ±2, ... is the 
mode number. 

The eigenvalue equation is easier to work with if re- 

1 /2 

expressed slightly. Solving for £,c{il) we find — 
2rj/{l + rf') = kx vac/^, and so 

/4:HN\ 

\avAc) 



V 



1 



where a — kx/k± < 1. If we now change variables to x 
defined by r/ = e"^ we find 



X 



ANH 



2Tm — i In 


TT 




tan — 


± TT 




7 





1 



7 



7 



coshx . 
(35) 

The variable x is convenient because physical quantities 
have simple expressions in terms of it. In particular, 
the phase velocity is Vph = (^i/k^ = wac Re coshx; 
the frequency's imaginary part is determined hy d = 
Im coshx/ Re coshx, and the position of the resonant 
layer is Zr = 2H In [Re cosh x] • 

Notice that Eq. i35\ takes the form x — ^ cosh x- For 
real A this equation has two real roots for x if ^ < Amax 
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Figure 2: Plot of the Wp/j /c^ — uji/kxCg against mode num- 
ber, n, where Cg is the adiabatic sound speed and uji is the 
real part of the eigenfrequencies, oj — uji{l + id), as pre- 
dicted by Eq. ( I35t . Different curves correspond to different 
background magnetic field strengths, Bq. Solid lines repre- 
sent resonances which are inside the radiative zone, whilst 
dashed lines correspond to unphysical modes whose reso- 
nances lie outside of the radiative zone. 



lOG 



and none if A > 



where Amax is defined as the 



maximum value taken by the function x/ coshx (which 
numerically occurs for x = Xmax = 1.1997). A similar 
result holds for A and x close to real, such as happens when 
uj = uji{l + id), and so solutions with d small only exist 
up to a maximum mode number, \n\ < n„iax, where 



2NH 



navAcJ coshxri 



Solutions with large d must be discarded because they lie 
beyond the scope of the linearized approximations within 
which we work. 

Figures 121 and |3] present our numerical solution of the 
eigenvalue spectrum obtained from Eq. ( I35t . Figure|2lplots 
Vph{n)/cs — u!i/{cskx) vs mode number n for various 
magnetic fields, Bq. Figure |3] similarly plots d against 
mode number for the same magnetic fields. In both plots 
the parameter a — Kx/K± is taken to be unity, a — 1. 
The result with a 7^ 1 is trivially found by rescaling Bq 
because Eq. i35\ depends only on the product aB^. In both 
figures a dashed line is used if the resonance of interest oc- 
curs above the top of the radiative zone (and so outside the 
domain of many of our approximations). 

Both of the figures |2l and |3l show two branches of solu- 
tions up to a maximum mode number, as expected. Their 
dependence on n can also be understood analytically if we 
take I X 1^ ~ 1, and so 77 ^ 1. In this limit we ob- 
tain from (I34> the following approximate equation for the 
MHD spectmm: 



UJi 

N 



2Ki 



nn 



{zr/2H + ln2) , 
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Figure 3: Plot of the d against mode number, n, where 
d gives the imaginary part of the eigenfrequencies, uj = 
0^1 (1 + id), as predicted by Eq. (I35t . Different curves cor- 
respond to different background magnetic field strengths, 
Bq. Solid lines represent resonances which are inside the 
radiative zone, whilst dashed lines correspond to unphys- 
ical modes whose resonances lie outside of the radiative 



d = 



In I tan(7r/7)| 



27rn(l - (zj2iJ + ln2)-i) ' 
which leads to the asymptotic solution 



wi 2K± , fmn 1 

— = In 

Al Trn \ avAc 

ln|tan(7r/7)| 
27rr7, 



(36) 



(37) 



These modes correspond to the branches of the figures |2] 
and|3]for which both Vph/cg and d fall with n. The other 
branch can be found by using the approximation |x| <C 1 
and hence cosh x « 1 + The spectrum in this case is 



UJi 



d = - 



aK±VAc 
H 
' avAc ' 



\2 


TT 


1 nn In 


tan — 




7 



(38) 



Note that for this branch Vph (n) / Cg is independent of the 
mode number, n, and d grows with n, as is also seen in the 
figure. 

In addition to requiring the resonance to occur in the ra- 
diative zone (the solid line in the figures), the validity of 
our approximations also demand the frequency to not be 
smaller than 10~^ s~^ due to our neglect of the 27 -day 
period solar rotation. Inspection of Fig.|2]shows that these 
two conditions are consistent with one another for a reason- 
ably large range of modes only for magnetic fields larger 
than a kG or so. 



4 The Resonant Properties 



Several striking features emerge from the previous sec- 
tion's analysis of the hydrodynamic eigenvalues and eigen- 
functions. The main two of these are: 

1. The eigenfrequencies predicted by Eq. i34\ are 
complex, implying both damped and exponentially- 
growing modes, and 

2. The eigenmodes are not smooth as functions of z 
about the singular resonant point, z — Zr{n, kx, ky), 
whose position depends on the quantum numbers of 
the mode in question. 

These features lie at the root of the surprisingly large 
effects which are produced by even moderate magnetic 
fields. We have argued that both of these features can be 
traced to a resonance between g-modes and Alfven waves, 
which arises due to the existence of gradients in the back- 
ground density, po{z)- 

In this section we expand on these two properties, to 
more precisely pin down their origin. In particular, we ex- 
pect the development of complex frequencies to indicate 
a hydrodynamic al instability. This is much Uke the imag- 
inary part which g-waves develop as one passes into the 
convective zone from the radiative zone. For g-modes this 
imaginary part arises because buoyancy no longer acts as a 
restoring force in the convective zone, but instead becomes 
destabilizing. Indeed the instability in this case drives the 
convection, which is why the radiative zone ends. We ex- 
pect a similar instability to arise for the Alfven/g-mode res- 
onance. 

4.1 Instability 

The necessity for complex frequencies is the most disturb- 
ing feature about the above eigenvalue problem, because it 
implies mode functions which grow exponentially in am- 
plitude with time. This signals an instability in the physics 
which pumps energy into these resonances, and this section 
aims to discuss the nature of this instability, and the physi- 
cal interpretation to be assigned to the imaginary part, d. 

Exponentially-growing instabilities within an approxi- 
mate (e.g. linearized) analysis reflect a system's propensity 
to leave the small-field regime, on which the validity of the 
approximate analysis relies. The question becomes: where 
does the instability lead, and what previously-neglected ef- 
fects stabilize the runaway behaviour? 

In the present instance we shall see that the normal 
modes become strongly peaked near the resonant radii, and 
that the energy flow near these radii is directed along the 
resonance plane, as would be true for a pure Alfven wave. 
This leads us to the following picture. Since helioseismic 
waves are likely generated by the turbulence at the bottom 
of the convective zone, it is natural to imagine starting the 
system with a regular helioseismic g-mode and asking how 
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it evolves. In this case the resonance allows the energy 
in this mode to be funnelled into the Alfven wave, and so 
to be channelled along the magnetic field lines away from 
the solar equatorial plane. This process drives us beyond 
the range of our approximations by making the wave am- 
plitude large near the resonant surface, but also by driving 
energy out of the equatorial plane where the rectangular 
analysis applies. 

We arrive in this way at an interpretation in which the 
instability signals the excitation of Alfven waves from g- 
modes due to their resonant mixing. In this picture the 
imaginary part of the frequency, F — Im u) = dui, de- 
scribes the rate at which the Alfven mode is excited in 
this process. Since the mixing and excitation are weak- 
field hydrodynamic phenomena it is reasonable that this 
rate can be computed using only the linearized magneto- 
hydrodynamical equations. 

Once excited, the mode amplitudes near the resonances 
grow until the energy in them becomes dissipated by ef- 
fects which are not captured by the approximate discussion 
we present here. The rate for this dissipation must grow as 
the mode amplitude grows, until it equals the production 
rate, F, at which point a steady state develops and the mode 
growth is stabilized. Because the resonant mode grows un- 
til its damping rate equals its production rate, it suffices 
to know the mode's growth rate in order to determine the 
overall on-resonance amplitude. 

How much more can be said about the amplitude of the 
modes depends on whether or not the mode stabilizes at 
small enough amplitudes to allow the hydrodynamic ap- 
proximation to the production rate, F = duii, to he ac- 
curate. For instance, if the mode only stabilizes once it 
is large enough to require a nonlinear analysis, then the 
final production and dissipation rates may be very differ- 
ent from the linearized expressions derived above. If, on 
the other hand, the modes saturate at comparatively small 
amplitudes by dissipating energy into non-hydrodynamical 
modes (such as by Landau damping), then it can happen 
that the stabilized mode amplitude is not large enough to 
significantly change the linearized prediction for its pro- 
duction rate. 

In what follows we use the imaginary part of oj to de- 
termine the resonant mode's amplitude at the position of 
the resonance. There is no loss in doing so to the extent 
that we take Im u; as a parameter which is not related to the 
mode quantum numbers by the hydrodynamical equations 
of motion. On occasion we shall also use the explicit for- 
mulae for Im OJ which is predicted by the hydrodynamics, 
and in this case we are assuming that the mode stabiliza- 
tion occurs at small enough amplitudes to not invalidate the 
assumption that this is a good approximation to the mode 
production rate. 

4.2 Spatial dependence 

We next address the singular shape of the wave fronts at the 
resonant points, z = Zr- Mathematically, this singularity 



arises because the interval of physical interest, < z < z*, 
corresponds to the interval < < ^ < and this con- 
tains the singular point ^ = 1 provided < 1 and ^* > 1. 
The hypergeometric functions can - and do - develop sin- 
gularities at ^ = 1, and it is this singular behaviour which 
we now examine. 

Using Eq. ( I67t of the Appendix B one finds as ^ ^ 1 
(from either side): 

62 oc 1^ — 11^" + nonsingular and -j^ oc |^ — l]""^, 

(39) 

with ly ~ (7 — l)/7 (giving 1^ = | for an ideal gas, for 
which 7 ~ |). For < < 1 we see that bz remains finite 
but dbz/dz ^ 00 as ^ 1. As we shall see, once the 
eigenfrequency's imaginary part is included db^/dz does 
not actually diverge as ^ — ^ 1, but instead adopts a system- 
atically large resonant line-shape. 

Using this information in eqs. (I14> allows the asymptotic 
form for all other quantities to be determined at the reso- 
nant point. In particular, since does not involve dbz/dz, 
it also remains finite on resonance, unlike all other com- 
ponents of the fluid velocity. This indicates that it is the 
transverse velocity perturbations (and so also the transverse 
kinetic energy of the oscillations) which increase most dra- 
matically at the resonant plane. 

Positions of resonances 

The equation for the position of resonances can be found 
by using the resonance condition 



Re ^ = Re £,c e 
and Eq. J36I 

Zs{n) = H 



z/H 



( kxVAc 



1 (40) 



\u!i{n) 

It takes the following transcendental form 



g..(„)/2H_in4 



(41) 

Taking into account the asymptotic form for the real part 
of the eigenfrequency, uji{n), as given by Eq. ( I37> . leads 
to the following asymptotic expression for the dependence 
of the resonant layer position on the node number, n: 



Zr{n) = 2H\ii 



here 



ANH 1 



Bn 



VAc 



In 



ANH 1 

cuVAc T^n 







cm s 



(42) 



(43) 



V43.4G^ 

is the Alfven velocity at the solar center 

Fig. 13 plots the position, Zr{n), predicted by Eq. ( I42> 
for the case of longitudinal wave propagation (ky = 0) 
and for magnetic fields in the range Bq =10G-lMGauss. 
(The same result for an obliquely-propagating wave with 
fcy ^ is produced by a larger value for Bq.) This plot is 
the analog of Figs. 5a,b in il4l . 

Spacing of resonances 
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Figure 4: A plot of the resonant position, Zr{n), vs 
mode number n for different magnetic fields in the range 




0.4 0.^ 
z/Rsun 



Figure 5: The distance between neighbouring Alfven res- 
onant layers vs the position of the layer within the solar 
interior 

Knowing the position of the resonances also permits us to 
determine the distance between them. This quantity is rel- 
evant to the propagation of particles like neutrinos through 
the resonant waves. The spacing is: 



+ ^^^^ 



1 



■ (44) 



This dependence of this quantity on mode number, n, is 
shown in Fig. |5]for longitudinal wave propagation, ky = 0, 
and for different magnetic field values. From this figure we 
see, in the region Zr > O.3i?0, that the distance between 
layer positions grows with magnetic field and with distance 
from the solar center. Both of these features were already 
seen in preliminary WKB calculations of the resonances 
(see Figs. 7a, 7b of ref. I.14J ). and our present numerics 
also confirm that the spacing is approximately proportional 
to _Bo for z > O.3i?0, that was seen earlier 

Eq. (I44> can be simplified using the expression, Eq. J42> . 



for the resonance position: 

Zr{n + 1) — Zr{n) 



2H 



(45) 



where the inverse mode number | n is proportional to 
the background magnetic field ^ Bq. This is in accordance 
with the spectrum Eq. ( I36t . for which the phase velocity 
Vpfi — LUi/kx equals the Alfven speed at the resonance 
position, uJi/tj,^ VA{zr,Bo) Bq. 

Numerically, it is noteworthy that the spacing between 
resonances can be in the ballpark of hundreds of kilome- 
ters. This is significant because this is close to the resonant 
oscillation length, for E ^ McV neutrinos, if - as 
now seems quite likely - resonant LMA oscillations are re- 
sponsible for explaining the solar neutrino problem since 



250 km (E/MeV) 
Ami sin 29 



(46) 



where Am| 



eV^. Repeatedly perturbing 
neutrinos over distance scales comparable to l^^^ has long 
been known to be a prerequisite for disturbing the stan- 
dard MS W picture of oscillations in the solar medium. This 
raises the possibility - recently explored in more detail in 
ref. 1 8 1 - that g-mode/ Alfven resonances can alter neutrino 
propagation. If so, the observation of resonant oscillations 
of solar neutrinos may provide some direct information 
about the properties of the MG waves we discuss here. 

Widths of the resonance layers 

We next examine the energy flow in the vicinity of a res- 
onance. This gives some information about the nature of 
the instability which is indicated by the complex mode fre- 
quencies. We find here that there is a large energy flow 
along the magnetic field lines, perpendicular to the z direc- 
tion. We use the dimension of the region in which this en- 
ergy flows as a measure of the width of the resonant layer 
In the linearized regime the energy flux, S, has both 
mechanical and electromagnetic (Poynting) components, 

S = Smech + Soni, with 



Smcch 



p V 

1 

An 



(47) 



{ [b' X [v X Bo]] } 



In order of magnitude, the relative size of these contri- 
butions is BqB' / {Airp') and this turns out to be small. 
This smallness is easy to see to the extent that B' /Bq ^ 
p' Ipo, since then BqB' / {inp') ^ Bq/{4:'kpq), which is 
small by virtue of the background inequality Bq / (Air) — 

Since the mechanical energy flux dominates, its direc- 
tion is simply given by the direction of v near the reso- 
nance. This may be determined using the asymptotic form 
as ^ ^ 1 in the expressions forp' and Vz in eqs. (I14t . No- 
tice that Vz approaches a finite value in this limit, while 
Vx, Vy and p' all diverge as ^ — > 1. This shows that the 
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energy flow, like the fluid flow, is mainly parallel to the 
resonant plane when evaluated near a resonance. 

Using the asymptotic expressions, dbz{z)/dz ~ S |1 — 
^|-'', and6^(z) |l-C|i-''+C, with constants C and 
1/ = (7- l)/7 « 0.4, we obtain cx |5^| ^B\l-i\^-'', 
and 



(48) 



p',Vx,Vy CX 






\B\ 


dz 






\B 






\B 




Cr2id\'' 





where we have used our definitions - uj = uji{l + id), 
with |d| < 1, and so ^ = ^^(1 - 2id), with = 

(jtxVAc/i^ij e^"^/^ - to write the real and imaginary parts 
explicitly. The energy flow near the resonant plane is sim- 
ilarly given by the time average S = (p'v) = iRe(p'v*), 
where the symbol (• • •) means average over time. Combin- 
ing the results for p' and v we see that near resonance 



1 



and so, for instance, Sz reaches half of its resonant value 
when ^ = ^± = l±dC where [(1 - + ^^±dY^^ = 
2{2dY. 

On the other hand, for z well beyond the resonant zone, 
where > 1, we know that bz{z) is given by Eq. i3Q\ . 
which states bz{z) oc ^-^^F{a,a*;l + 2Kj_;^^'^)- By 
virtue of the condition a* = 1 — c + a we see that this 
is a real function (in the limit d ^ 0). Consequently, 
the explicit factors of 'i' in eqs. il4\ imply that quanti- 
ties like the pressure perturbation, p'{z), and the trans- 
verse velocities Vx,y{z) are pure imaginary. It therefore 
follows that the z-component of the energy flux vanishes 
in the region z > Zr, outside of the resonance, because 
5*2 — ^Ke{p'v*) — 0, while the transverse components 
are nonzero in this region, but are falling off exponentially. 

This peaking, then vanishing, of Sz motivates our choice 
of the width of the Sz{z) resonant line-shape as a measure 
of the resonance layer width. From this condition we can 
easy find the results for the widths z± (or, equivalently, 

(50) 



giving our final expression for the width of the Alfven sin- 
gular layer: 

6z = z+-z^= Hln (^^^ w 4iJ|d|2'^/('^-i) . (51) 

Substituting the expression for d in terms of mode number, 
using Eq. ( I37t . this becomes 



Azrin) w — 2'^/('^-i) In 
nn 



tan ■ 



(52) 
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Figure 6: A plot of neighboring density profiles, show- 
ing their width in comparison with the spacing in be- 
tween them. The plot is done for resonances in the region 

Zr ^ O.3i?0. 

Together with Eq. ( 145 > we then get the ratio of the dis- 
tance between resonances to the width of each resonance. 



Zr{n + 1) — Zr{n) 

Azr{n) 



27/(7-1) In I tan ^ I 



0.5. (53) 



Figure|6]plots the density profiles of two neighboring reso- 
nances in the region Zr ^ 0.3Rq. In this region their width 
is approximately twice their separation as measured by the 
distance between peaks. That is not surprising since the 
shape of the density profile is the same as in Eq. (I48t for 
other MHD perturbations (see below, Eq. (I55» . 

Density Profiles 

We now give the explicit expression for the density pro- 
file predicted for our wave solutions. The basic expression 
is obtained through the direct substitution of our expres- 
sion for bz{z) - keeping in mind eqs. J24> . ( I28> and i29\ 
- into the last of eqs. ( I14> . We find in this way the fol- 
lowing expression for p' / po, divided by the arbitrary wave 
normalization, Di: 



Di \po{z) 



7 



klcl 



7 

.2 _ 1^2 



1,2 2 



a2,^2;c2;C) 



F{ai,bi\ci;C) 

F{ai,bi;ci;£,c) 



F(a2, 62; C2; fc) 



1 - ip V6 

F(ai,6i;ci;^c 



i^(ai + l,6i + l;ci + l;0 



F(a2 + l,62 + l;c2 + l;e) 



F(a2,62;c2;^c) 



(54) 



where ai is given by Eq. M2\ for q — ~i(3, (T2 = (7 — 
l)/27 - i/3/2; 01,2 = en, 2 + -ft'j., &i,2 = cri,2 - 
ci,2 = l±i(3 (see Eq. (|29}). 



13 



The behaviour of the density perturbation near the reso- 
nant point, Zr, takes a simpler form for the Lagrange per- 
turbation 

5p{f) — p'{f) + 5r ■ \/po = 

= [p'{z) - {Sz/H)po{z)]e-'^^'-''^--''yy\ (55) 

where the Eulerian perturbation p'{z) is given by Eq. ( I54> . 
In linear MHD the Lagrangian density perturbation 
is directly constrained by mass conservation, Eq. Q, 
—iujp^^Sp{z) — — u, where u = divv is the compress- 
ibility, given by Eq. ( I15I I. One obtains in this way the fol- 
lowing Lagrange density perturbation near the resonance, 
z ^ Zr{n): 



1 Spiz) 



Di po 
where 



-(7-l)/7 



(56) 



C = 



47 



7- 



Here the phase velocity, Vph{n) = ijji{n)/kx, is given by 
the solution to the eigen-spectrum, oj = a;i(l + id), such 
as obtained in Eq. (( I37» . and we include the unknown nor- 
malization factor Di. 

This quantity \5p/ pq] is plotted in Fig.|6lfor two neigh- 
bouring resonances near z — O.Si?© for the two values 
of the magnetic field, Bq =10kG and Bq =100kG. We 
see that the presence of the resonances introduces a se- 
ries of density excursions, located at the resonant posi- 
tions. The disappearance of the resonant density spikes 
in the zero-field limit is most easily seen by noticing that 
in the limit i?o ^ the position z^ becomes negative, 
Zr = —H\n4 < (see Eq. 1411}. or occurs outside the 
physical region < z < z*. 

5 Possibilities for Observations 

Given the strength of the resonant effects on the heUoseis- 
mic wave profiles, we next explore some of their observ- 
able properties to determine how they might be detected, 
and to what extent their effects are consistent with current 
heUoseismic observations. 

5.1 Measurements of Mode Frequencies? 

Many helioseismic waves have been detected, but to date 
all of these are pressure-driven p-modes rather than the 
usually-lower-frequency g-modes. The latter are more dif- 
ficult to detect because they are so strongly attenuated 
within the convective zone. 

Because the central Alfven frequency is so small com- 
pared with the other helioseismic frequencies, we only ex- 
pect resonances between Alfven waves and those helioseis- 
mic waves which penetrate most deeply into the solar in- 
terior, and which have the lowest frequencies. This is why 



the resonance we have found involves mixing only with g- 
modes, and not p-modes. It follows that so long as observa- 
tions are restricted to p-modes, the resonances we describe 
here will have negligible implications for the comparison 
between measured p-wave frequencies and solar models. 

Things become more interesting should g-modes be ob- 
served, since then a comparison between observations and 
predicted frequencies could depend more strongly on the 
solar magnetic field. Figure (|2j (or Eqs. ( I37> and i42l ) 
show how the frequencies of the resonant modes depend 
on magnetic field. Unfortunately, frequencies this low are 
likely to be a challenge to detect for quite some time. 

Of course, a proper comparison requires the inclusion 
of the full (spherical) geometry of the problem, as well as 
temperature and gravity gradients, solar rotation, and other 
factors which we ignore here. Our goal in this paper is 
merely to point out the potential for larger-than-expected 
effects, and to obtain a rough first determination of the sen- 
sitivity to magnetic fields. 

5.2 Constraints on Mode Amplitudes 

Mode amplitude remains undetermined within the linear 
approximations assumed here, and is ultimately controlled 
by the mechanisms which stimulate and damp the waves of 
interest in the solar environment. Amplitude can neverthe- 
less be characterized in any of three equivalent ways: the 
size of the mode's contribution to the mass-density profile, 
the size of the fluid velocity, v^, at the solar surface, and 
the total mode energy. In principle, our solutions can be 
used to express any two of these in terms of the third, as 
functions of the background magnetic field. 

A sufficiently detailed comparison of the resulting ex- 
pressions - should this become possible - would then al- 
low one to obtain the magnetic field strength. We illustrate 
this by calculating here the energy content of the resonant 
MG waves. 

Mode Energy 

The mode energy is a useful way to quantify the mode am- 
plitude. Here we provide formulae which make this con- 
nection explicit. The mode energy is most simply com- 
puted by using the virial theorem, which relates the total 
energy to the kinetic energy averaged over a single period 
of oscillations: E — 2(i?kin): 

(£^ki„) = ^ / Po ([Re v(x, t)f) . (57) 

For modes having a small imaginary frequency, \d\ ^ 1, 
one obtains the following simple approximate formula for 
the above average 



([Rev(x,i)]2) = i|v(z)| 



(58) 



+ -[(lmt;,(z))2-(Re«,(z))2] 
X sm[2{kxX + kyu)] + 0{d) , 
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The second term in this expression is typically negligible 
once it is integrated over dx and dy. 

Using our low frequency approximation, lu^ <C fc^c^, 
we find the following result for the squared velocity. 



12. 



MO 



Kim- 



d^ 



dC 



(59) 



We next integrate this result over the solar volume. In 
our simplified geometry we do so by bounding the trans- 
verse directions with a circle of radius Rq, having cross- 
sectional area ttRq. This gives the basic expression for the 
energy per unit length as dE/dz = (7r/4)i?Qpo(2:)|v(z)p. 

A useful expression is obtained by eliminating the un- 
known normalization, Di between eqs. ( 15 9> and (I54> . 
This gives a direct, but complicated, relation between the 
mode energy and the size of the peak density excursion 
on resonance. To this end we denote the function in 
brackets in Eq. ( I59> by v4(f). With this choice we have 
|v|' = Vphin)A{t n), with vl^{n) = co^kl = v\{zr) = 
j^2^Zr{n)/H I (47rp^). Note that this last expression uses the 
resonance condition, ^ = {kxyAizr))^ /^"^ = 1, for the 
mode eigenfrequencies. If we similarly denote the right- 
hand-side of Eq. (|53 by B{Cj, so p' / po ^ Di then 
we find 



E 



MHD 



(n) = So 



-V 



1 



'Re p'{zr{n)) 
po{zr{n)) 



1 Gauss/ (ReB(C = l,n))2 
' '■^'d^A{^,n) 



e Dt{n) 



(60) 



where for convenience we assume Di to be real. In this 
expression we have used the explicit background density 
distribution, po — Pc e~^/^, and changed variables from 
z to ^. The lower limit of integration is — e^^-^'^'^^l^ 
which takes values 1 > ^ 10^"^ for a resonance between 
the solar center and the bottom of the convective zone. We 
take as upper limit ^* = e("'^^0^^' )/^ 3> 1 for the same 
choice of the Alfven resonance positions. 

The energy scale, Eq, in Eq. ( I60t denotes the combina- 
tion Eq = [1/1&)R%H X (1 Gauss)2. Using (1 Gauss)^ = 
1 erg/cm^'^ gives Eq ^ 2.1 x 10"^" erg. From this we can 
see that, for magnetic fields Bq — 10^ — 10'' Gauss, density 
spikes as large as p' / po ~ 1 — 10 per cent will involve ener- 
gies E < 10'^^ erg, provided A/ and B are order unity, 
which is reasonable for the maximum g-mode energy. 

Density Perturbations and Neutrino Oscillations 

The g-mode/Alfven resonances could give rise density per- 
turbations with a large amplitude. Are such density excur- 
sions constrained by observations? At first sight one might 
think so, because the success of the comparison of standard 
solar models with helioseismic observations constrains the 



solar density profile to be within 1 per cent of solar model 
predictions. Since these models ignore magnetic fields, one 
might expect magnetic-field-induced effects like those we 
consider here to be ruled out at this level. 

Unfortunately, deconvoluting solar properties from ob- 
served wave profiles is at present only possible subject to 
smoothness assumptions. As a result current constraints on 
deviations of p from solar model predictions only apply to 
profiles which do not vary on scales shorter than around 
1000 km (9l[T0|[ni. It is for this reason that present anal- 
yses do not yet exclude the existence of such resonances at 
an interesting level. Observable effects are only possible 
if the distance between spikes is comparable to the neu- 
trino oscillation length at the position where neutrino con- 
versions occur Furthermore, their amplitude must be of 
order a few percent in order to produce observable effects. 

If these density spikes are sufficiently large, they could 
have implications for neutrino oscillations, and we have 
argued 1 8 1 that this may be the first place where they are 
potentially detectable. Besides identifying whether the re- 
quired density profiles can be obtained using reasonable 
magnetic fields, a proper analysis must also ask whether 
the resulting modification of the signal can be identified 
within the considerable uncertainties which are inherent in 
any neutrino measurement. 

The sensitivity of solar neutrino oscillations to noise in 
the solar interior has been re-examined in ref. |8|, using 
the best current estimates of neutrino properties. The out- 
come is that the measurement of neutrino properties at 
KamLAND now provides potentially important informa- 
tion on fluctuations in the solar environment on scales VJ^^ 
(see Eq. ( !46» to which standard heUoseismic constraints 
are currently insensitive. 

Performing a combined fit of KamLAND and solar data 
it has also been shown how the determination of neutrino 
oscillation parameters strongly depends on the magnitude 
of solar density fluctuations. 

The Alfven/y-mode resonances turn out to have the right 
spacing and position for magnetic fields in the range of 
So = 10 — 100 kG. Depending on the size of the magnetic 
fields which are ultimately found in the radiative zone, this 
perhaps opens a completely new and interesting observa- 
tional window on the solar interior. 



6 Discussion 

Here we have presented an analytic discussion of the influ- 
ence of magnetic fields on heUoseismic waves in a simpli- 
fied geometry. We have argued that this geometry captures 
the physics not too close to the solar center (r > H ^ 
O.li?0) within the radiative zone, provided that the back- 
ground magnetic field is perpendicular to the background 
density gradient. In particular it would apply to a dipole 
magnetic field configuration near the dipole's equatorial 
plane. 

We find that within our approximations sufficiently large 
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magnetic fields can cause significant changes to the profiles 
of helioseismic g-waves, while not appreciably perturbing 
helioseismic p-waves. The comparatively large 5-wave ef- 
fects arise because of a resonance which occurs between 
the g-modes and magnetic Alfven waves in the solar ra- 
diative zone. Although the radiative-zone magnetic fields 
required to produce observable effects are larger than are 
often considered - more than a few kG - they are not di- 
rectly ruled out by any observations. Noting and character- 
izing the basic features of such resonance constitutes our 
main result. 

Although the density profiles at their maxima could have 
amplitudes as large as a few percent or more on resonance, 
we do not believe the corresponding radiative-zone mag- 
netic fields can yet be ruled out by comparison with helio- 
seismic data, since the density excursions are sufficiently 
narrow (hundreds of kilometers) as to evade the assump- 
tions which underlie standard helioseismic analyses. 

For magnetic fields in the 10 kG range, the best hopes for 
detection of the resonant waves may be through their influ- 
ence on neutrino propagation. This influence essentially 
arises because the presence of strong density variations af- 
fects the solar neutrino survival probability, with a corre- 
sponding change in the resulting solar neutrino fluxes. As 
shown in ref. lHJ, the measurement of neutrino properties 
at KamLAND provides new information about fluctuations 
in the solar environment on correlation length scales close 
to 100 km, to which standard helioseismic constraints are 
largely insensitive. It has also been shown how the determi- 
nation of neutrino oscillation parameters from a combined 
fit of KamLAND and solar data depends strongly on the 
magnitude of solar density fluctuations. 

Since the resonances rely on the condition that B _L V/o, 
there are several magnetic-field geometries to which our 
analysis might apply, of which we consider here two il- 
lustrative extremes. Suppose first that, in spherical coordi- 
nates (r, 9, (p), we imagine Br ^ but Be 7^ 0. Then the 
field is always perpendicular to a radial density gradient 
and the resonance we find might be expected to arise in all 
directions as one comes away from the solar center In this 
case the solar g-modes would tend to be trapped behind 
the resonance, and so are kept away from the solar surface 
even more strongly than is normally expected. This would 
make the prospects for their eventual detection very poor. 

Alternatively, if the magnetic field has more of a dipole 
form it might be imagined that the condition B _L Vp only 
holds near the solar equatorial plane, and not near the solar 
poles. In this case a more detailed calculation is needed in 
order to determine the resulting wave form. This kind of 
geometry could have interesting consequences for the so- 
lar neutrino signal, because in this case the deviation from 
standard MSW analyses only arises for neutrinos which 
travel sufficiently close to the solar equatorial plane. Given 
the roughly 7-degree inclination of the Earth's orbit relative 
to the plane of the solar equator, there is a possibility of ob- 
serving a seasonal dependence in the observed solar neu- 
trino flux. Since the presence of the MG resonance tends 



to decrease the MSW effect, the prediction would be that 
the observed rate of solar electron-neutrino events is maxi- 
mized when the Earth is closest to the solar equatorial plane 
(December and June) and is minimized when furthest from 
this plane (March and September). 
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Appendix A. Some Hypergeometric 
Facts 

We here record some properties of Hypergeometric func- 
tions which are useful in obtaining our results. The series 
solution to the Hypergeometric equation, Eq. i2Q\ . which 
is analytic throughout the unit disk |^| < 1 is: 



F{a,b;c;0 - 1- 



ab 



a{a+l)b{b+l) 
c(c+l) 



(61) 

for c 7^ 0,-1,-2,.... The function is defined elsewhere 
in the complex ^ plane by analytic continuation from this 
series expression. 

A convenient way to describe the properties of F for 
1^1 > 1 is to interchange the interior and exterior of the unit 
circle using the mapping ^ ^ 1/^. Because this is a par- 
ticular element of the Hypergeometric group, which maps 
the three singular points of Eq. i2Q\ (^ — 0, 1, 00) among 
themselves, it has a simple representation on the functions 
F(a, b; c; ^), given explicitly by the following identity 0: 



F{a,b;c;0 
r(c) 
r{b-a) 



T{b)T{c-a) 
r(a - b) 

'r{a)r{c-b) 



(62) 

F{a,l-c + a;l-b + a; C^) + 
(-C)-^ F (6, 1 - c + &; 1 - a + 6; T') • 



This identity often can be usefully combined with the series 
expression, Eq. i61\ to obtain asymptotic forms as |^| — > 

00. 

The series solutions given by eqs. ( 16 1> . i62\ are calcu- 
lated numerically very slowly from both sides in the vicin- 
ity of the Alfven resonance, ^ < 1 and ^ > 1 respectively. 
For these regions we have used other representations of the 
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Hypergeometric functions For^ < 1: 

F{a,b;c;£_) r{c-a-b) 



r(c) T{c- a)T{c-b) 

F {a, b;a + b - c + l;l - £_) 

+ r(a)r(5) 

xi^(c — a, c — 6;c— a-6+l;l — ^) 



Appendix B. Asymptotic Forms 

The identities of the preceding Appendix A can be used to 
derive the asymptotic forms in the hmit of large ^ for the 
solutions Yi and Y2 discussed in the main text. For these 
purposes we use the following more convenient basis of 
solutions to Eq. (I2Q> : 



(63) 



and for ^ > 1: 

F(a,6;c;0 



Y3 = C'^F (6+l-c,6;a + 6- c+l;l-i 



T{c-a-b) 



z—a — b 



r(c) T{c-a)T{c-b) 
F(a,a- c + l;a + 6 - c + 1;1 - C^) 
r(a + 6-c) 

+ r(a)r(6) ^ 

F (c — a,l — a; c — a — 



y.F [ c — aA ~ a;c ~ a 



(67) 



— a — b \^ 



l;l-r') ■ (64) 



Another asymptotic form of the Hypergeometric func- 
tion which proves useful for our purposes is the second 
Watson's form for large values of the parameters a, b and 
finite c. We use this in its following variant (Sj: 

F{a + \,b- X;c;^- |) = 



The expansion of our original solutions, Yi and Y2, in 
terms of these new ones has the form Yi = NiY^ + N2Y4 
and Y2 ^ N[Y3 + I4, with coefficients iVi,2 and iV( 2 
given by 

r(c)r(c-a-6) 



N2 - 



r(i -6 + A)r(c) 



r{c-a)T{c-b) ' 
T{c)r{a + b- c) 

r(a)r(6) ' 

r{2- c)T{c - a- b) 

r(i-a)r(i-6) 



r(i/2)r(c-6 

xii + uy-"-'' 



A) 

5 X A-1/2 



2 X 



r(2-c)r(a + b-c) 

= r(a + i-c)r(6 + i-c) • ^^^^ 

The large-^ limit is now easily obtained using the result 



x{l + Oi\X\-')) 



(65) 



F(a,6;c;l) 



r(c)r(c-a-6) 



Here the argument is connected with the auxiliary vari- 
able z = 1 — 2^ and corresponds to u = [z + Vz^ — 
(if Im z < 0) and to u = [z + \J z^ - 1] (if Im z > 0) 
with the same rule for the choice of upper (lower) sign in 
the exponential factor within brackets. 

The correction 0(A^^) <C 1 turns out to be important 
for the case A = Ky^ ^ 1 which can be used for high 
frequency spectrum oji < N (although we do not present 
this result here). We find the following correction in IrSTI : 

2A 1 _|_ g±i7r(c-i)^a-()+2A ' 

1 L + Mu-^ + Nu-^ 

2 1 ' 
1L + MU + Nv? 

" 2 ' 

L={a + b-2c+lf-a + b-^, 

M = -2{a + 6 - l)(a + 6 - 2c + 1) , 

1 



T{c- a)T{c-b) ' 

which is finite provided Re {c — a ~ b) > and c 7^ 
0, —1, —2, ... |5|. Notice that these conditions are satisfied 
for the choices of a. 6 and c obtained in the main text. Ap- 
plying these expressions to Y3 and Y4 then gives, in the 
limit ^ — + CX3: 



Y. 



Y4 



_^T{l + a + b-c)r{a~b) 



r(a)r(a + 1 -c) 



_,_,r{l- a- b + c)T{a-b) 



(69) 



r(i -5)r(c-5) 

Using the previously-given expression, Eq. ( I68t . for the co- 
efficients a, b and c, the limiting form for original solutions 
^1.2(0 finally become: 

_,r(c)r(a-6)^ 



r(a)r(c- 6) 



Yo 



r(2-c)r(a-b) , 



N={a + b-lY + a-b+- 



(66) 



r(i -6)r(n-a-c) 

with coefficients 2 given by the elementary functions 



(70) 



We use this expression obtaining the low frequency 
spectrum J34> . uji{n) ^ A^, and elsewhere in the deriva- 
tion of approximate forms. 



h 



sin7r(c — a) — (— ) 



sin7r(c —a — b) 



simra-i-Y " ^sim:{c-b) 

h = : 1 -Ts • (71) 

sin7r(c —a — b) 
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Notice that these quantities satisfy the relation I2/I1 = 

_„iTrc 
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